Pairing of Solitons in Two-Dimensional S = 1 Magnets 
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We discuss the structure of topological solitons in a general non-Heisenberg model of isotropic two- 
dimensional magnet with spin 5* = 1, in the vicinity of a special point where the model symmetry 
is enhanced to 5*17(3). It is shown that upon perturbing the 5'L*'(3) symmetry, solitons with odd 
topological charge become unstable and bind into pairs. 
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Introduction. - Many condensed matter systems can be 
successfully described with the help of effective contin- 
uum field models. In systems with reduced spatial di- 
mensionality, topologically nontrivial field configurations 
are known to play an important role Magnetic sys- 
tems are usually modeled with the help of the Heisenberg 
exchange interaction. In many instances the fluctuations 
of the length of the local magnetic moment occur at a 
large energy scale and can be neglected; the continuum 
field theory in that case is the so-called 0(3) nonlinear 
sigma model (NLSM) describing the dynamics of a three- 
component real unit vector field, and the topological ex- 
citations of this model are well understood 

However, for spin S > 1/2 the general isotropic ex- 
change goes beyond the purely Heisenberg interaction bi- 
linear in spin operators Si , and may include higher-order 
terms of the type (SiSj)^ with n up to 25*. Particularly, a 
general 5* = 1 model with the isotropic nearest-neighbor 
exchange on a two-dimensional (2d) square lattice is de- 
scribed by the Hamiltonian 



write down the spin-1 state at site j as 



n 



J(S,S,)+if (S,S, 



where (ij) denotes the sum over nearest neighbors, and 
J and K > are respectively the bilinear (Heisenberg) 
and biquadratic exchange constants. The model U]) has 
been discussed recently in connection with 5=1 bosonic 
gases in optical lattices Q and in the context of the de- 
confined quantum criticality 0, [1]. The effective field 
theory for the above model is generally more complicated 
than NLSM: the order parameter belongs to the 2d com- 
plex projective space CP^, and at two points, J = K 
and J = 0, the model symmetry is enlarged to SU{3). 
The aim of the present paper is to show that the crossover 
from SU (3) to SU{2) symmetry, taking place in the vicin- 
ity of those special points, features drastic changes in the 
structure of topological excitations, which correspond to 
pairing of solitons of the CP^ model. 

Continuum field description.- The spin-1 state \'4')j at 
a given site j is a linear superposition of three basis states 

cr = 0,±1. It is convenient to 



(2) 



using the "cartesian" states \z) = |0), |x) = (| — 1) — | + 
l))/V2, |y) = i(| - 1) + I + l))/\/2, then the three num- 
bers tja transform under rotations as the components of 
a complex vector tj . The normalization condition brings 
the constraint 't*j - tj — 1 • The states ([2]) can be viewed as 
^[/(B) coherent states corresponding to the bosonic op- 
erators ij,a, and the S = \ operator can be represented 
as 5° — —i^ahc^j bij c- Taking into account that the state 
(PI) may corrtain an arbitrary overall phase factor, one 
concludes that the order parameter space M of the prob- 
lem is four-dimensional and isomorphic to CP^. 

The lattice Lagrangian of the model expressed in terms 
of the complex unit vector t takes the form 



(3) 



(1) where the local Hamiltonian average (/li.j) is given by 



{\,) = j{t: ■ t^){t] ■ t,) + ( J - K)it* ■ t*)it, ■ t^). (4) 

This makes obvious that the system is always invariant 
under global rotations tj.a ^ 'Ti-abtj,b, with an arbitrary 
0(3) rotation matrix TZ, as well as under local "gauge" 
transformation tj i— > tje^^^ . At J = K the symmetry 
becomes higher as there is an invariance under a global 
transformation tj i-^ Utj, with U G SU{3). Moreover, if 
the lattice is bipartite, at J = the energy is invariant 
under making an arbitrary SU (3) rotation on the sites 
belonging to one sublattice if this is accompanied by a 
conjugate transformation tj i-^ U*tj at the other sublat- 
tice, so the point J = is iS'J7(3)-invariant as well. 

Breaking up the complex vector t = u + iv into two 
real vectors representing its real and imaginary parts, one 
can write the on-site spin and quadrupole averages as 



\a)j with Sj\a)j — cr\a 
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(S) = 2{u X v), 

Sab = {SaSb + SbSa) 



2((5ab 



UaUb - VaVb 



(5) 
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One can use a different parametrization, directly con- 
nected to the physical averages, by introducing the eight- 
component vector n, 



tr(t* • Xat), 



(6) 



where Aq, a = 1, . . . 8 are the well-known Gell-Mann ma- 
trices which form, together with a unit matrix 1, a basis 
in the SU{3) matrix space. The vector n is subject to 
the following two constraints: 



n2=4/3, n • (n*n) =8/(3V3), 



(7) 



where the *-product of any two vectors n and n' is de- 
fined as {n*n')a = V^dafd-yn^n'^, and dajS-y are the struc- 
ture constants defined by the anticommutation properties 
of the Gell-Mann matrices, {Aq, X/j} = ^Sa/3t + 2daf3'yX'y 
One can show that the constraints ([7]) in fact reduce the 
dimension of the n-space to four. The quantities Ua cor- 
respond to the following on-site averages: 

714 = Sxz, = Syz, rii = Sxy (8) 

n^ = {Sx.- Syy)/2, n8 = \/3(5,j2-2/3), 

which can be split into the vector of spin averages m and 
the vector of quadrupolar averages d, 

m = (ny, -715,712), d = (711,713, 714,716, 7i8). (9) 
In those variables, the Hamiltonian takes the simple form 

{Kj) = - • d,) + - 2J){m, ■ m,), (10) 

which explicitly shows that J > K corresponds to a fer- 
romagnet (FM), J < to an antifcrromagnct (AFM), 
and < J < -ftT to a quadrupolar (spin nematic) order 
(hereafter we assume that K > and will not discuss the 
so-called orthogonal spin nematic present at A' < 0). 

In terms of n, the lattice Lagrangian can be written as 
C = J2j ^{''T'j) ~ (^y), with the dynamic part 



^, , 3no • {nAdtn) 
$(n) = -- 
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(11) 



Here the 5'C/(3)-crossproduct is defined as (n A n')^ = 
fafS-ynfju'^, where fa/s^y is another set of structure con- 
stants defined by commutators of the group generators 
[Xa, A/3] — 2ifa/3'yX-y, and no is an arbitrary vector satis- 
fying the constraints ([7]). 

Topological analysis. - To describe topological solitons, 
one needs to pass to the continuum description first. The 
continuum Lagrangian of the model ^ can be obtained 
by the gradient expansion of the discrete energy W re- 
taining the leading terms, that gives W = J d^xw with 

w^J{\d^t\^-\t*-d^t\^} + {J-K)\t^\^ 

-{J-K)[\t- d,t\^ + \ [t'id.tf + c.c] }, (12) 



where /z runs over space coordinates {x,y). The above 
form is valid for the region J > K/2, where the short- 
range spin-spin correlations are of the ferromagnetic 
type, as can be seen from pHI) . 

To classify the topological excitations, one needs to 
know the so-called degeneracy space M^i that includes 
all values of the order parameter field corresponding to 
the ground state of the system. For the model ^ the 
space M/j is continuous and depends on the type of the 
ground state: for FM or AFM it coincides with the unit 
sphere 5^, for the nematic case it is a 2d real projective 
space RP^ = S'^/Z'^ (a unit sphere with the opposite 
points identified), and at J = K the degeneracy space 
is enlarged to CP^. For all the above spaces, the sec- 
ond homotopy group is nontrivial, 7r2(M£)) = Z which 
makes possible the existence of so-called localized topo- 
logical solitons, whose order parameter distribution be- 
comes uniform away from some point. 

If the order parameter lies completely in M^i, the en- 
ergy contains only terms with gradients, so there is no 
natural space scale. If corresponding soliton solutions ex- 
ist, they have a finite energy which does not depend on 
their size, and are stable against collapse. Another possi- 
bility is to allow the order parameter to leave M^, which 
breaks the scale invariance. Static solitons of that type 
are unstable against collapse due to the Hobart-Derrick 
theorem, but they can be stabilized by some internal dy- 
namics P, Q • We will study the structure of both types 
of solitons for the model ([1]). 

For the sake of analyzing static soliton solutions the 
Lagrangian ([3]) with the energy (fT2|) is equivalent to the 
2d CP^ model 6] with an additional "anisotropy term" 
proportional to (J — K). Let us start from the SU{i)- 
symmetric point J = K. In that case a localized topo- 
logical soliton corresponds to the field configuration with 
nonzero topological charge Q: 



2-K 



d'^xe^^{d^t* ■ dyt), 



(13) 



where the indices ^, v run over {x,y). The invariant lie 
takes only integer values and corresponds to the mapping 
of the compactified 2d space onto CP^. The exact 
q ~ 1 soliton solution is well known [6|: 



t={^a + zb)/^/WTe 



(14) 



where z ~ x + iy is the complex coordinate (the soliton 
center is assumed to be at the origin), a and b are two 
mutually orthonormal complex vectors, and ^ has the 
meaning of the soliton size. The energy of such excitation 
according to (fT2)) is E — 2ttK . For an arbitrary value of 
q, the general soliton solution can be written as 



ta = 



la 



(Eal/aP)^/^ 



2^/2 J /a — Ca (z Zk,a)j 



k=l 



a = x,y,z, 
(15) 
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and the corresponding energy is E = 2Tr K\q\. 

Ferromagnetic solitons.- On the ferromagnetic side 
J > K the minimum of energy is achieved for 



t = (ei +ie2)/V2 



(16) 



with ei_2 being a pair of orthogonal real unit vectors. In 
that case on the degeneracy space Md the order param- 
eter is equivalent to the unit vector m = (ei x 62) (a 
rotation around m corresponds to a change of the over- 
all phase factor t 1-^ te*''' and thus does not change the 
physical state). Thus, localized topological solitons for 
J > K correspond to the mapping i— > and are 
characterized by another topological charge 



(17) 



It is easy to calculate the topological charge for a re- 
stricted field configuration satisfying (jl6p : a general pair 
of orthonormal vectors ei.2 can be obtained from e^^y by 
an arbitrary rotation TZ(0, ip^-tp), where and ip are re- 
spectively the polar and azimuthal angles characterizing 
the direction of the unit magnetization vector m, and 
the third angle ip corresponds to the rotation around m. 
A straightforward calculation yields 



1 

2^ 



(18) 



One is led to conclude that solitons of the CP^ model 
tend to pair upon perturbing the SU{3) symmetry, which 
constitutes the central observation of the present paper. 

The above result can be also obtained by noticing 
that for the configurations the energy takes the 

form W = (J/2) / (Px {d^m)'^ . This is exactly the en- 
ergy of the 0(3) NLSM, and the well-known Belavin- 
Polyakov (BP) soliton solution with the topological 
charge Qm — 1 will have the energy E =^ At: J, which in 
the limit J K is twice the energy of the q — 1 soliton 
of the CP^ model. In fact, one can explicitly check 
that the ferromagnetic BP soliton is a particular case of 
the general solution (fTSjl with q — 2. 

Solitons for spin nematic. - On the nematic side J < K 
the minimum of energy is reached for t = ue^^ , where u is 
a real unit vector and x is an arbitrary phase. The degen- 
eracy space is thus = RP^. The energy then takes 
the form W = K J cPx (9^it)^, where u must be under- 
stood as a director, i.e., u and —u are physically identi- 
cal. It is worth noting that in contrast to the other phases 
the spin nematic allows for a nontrivial 7ri-topological 
charge as well, ni{RP'^) — Z2. 

If one defines the topological charge Qu according to 
p?)) . simply replacing m by u, then in the BP soliton 
with Qu — 1 the director u goes over Md twice; the 
energy of such a solution is Ebp — St: J. However, the 
director property of u allows one to construct a solution 
with u going over Md just once, which has Qu = h 



and the energy Ebp = An J. In the limit J ^ K this is 
again twice as much as the energy of the q = 1 solution 
Pl)) . which suggests that this soliton is a descendant of 
the q = 2 solution of the CP^ model. This indicates that 
the tendency to pairing exists on the nematic side as well. 

The fate of solitons with q ~ l.~ Up to now we have 
considered only static solitons with the order parameter 
lying completely inside M^. We found that ioi J—K 
the lowest energy solutions of that type are descendants 
oi q — 2 soliton of the CP^ model, while the q = 1 solu- 
tion seems to exist only at J = K. To get further under- 
standing of what happens in the vicinity of the SU (3)- 
symmetric point J = K, let us discuss the CP^-soliton 
with g = 1 for small but finite J — K. One can easily see 
that at J — K ^ any solutions with q = 1 must involve 
a deviation of the order parameter from the degeneracy 
space Md. Due to the Hobbart-Derrick theorem, this 
means instability of static solitons with q — 1 against 
collapse. However, q = 1 solitons can be stabilized by 
internal dynamics in presence of additional integrals of 
motion, e.g., stable solitons with the magnetization vec- 
tor precessingaround the easy axis exist in the uniaxial 
ferromagnet In our case, it is also possible to con- 
struct such a solution. In terms of the complex vector 
t ^ u + iv this is a planar configuration, where u and 
V are parallel to the plane (1, 2) orthogonal to some axis 
63, for definiteness let it be the z axis (a more general 
solution can be obtained by an arbitrary rotation). It 
is convenient to use the 8- vector notation ([6]): only four 
components of n are nonzero and it takes the form 



n= {R^,R,,Ry, 0,0, 0,0, 1/V3), 



(19) 



where i2 is a unit vector combining one spin average R^ = 
TO3 and two quadrupolar variables R^ = di, Ry = d2 
(cf. Using PH)) and pT|) . one obtains the efi^ective 

Lagrangian for the chosen subspace. 



Wr 



1 \ - Rq ■ 

2 ^ 1 



Rq ■ {Rj X dtRj) 



Rq • Rj 



Wr 



(20) 



^ [ — R^RJ 

<ij> 



{J ~ K)R,^,R,^j\ 



where Rq = (0,0,-1), and in pT|) we have used no — 
(0,0,-1,0,0,0,0,^). The Lagrangian ^ describes 
the dynamics of a classical anisotropic ferromagnet with 
the unit magnetization vector R; the anisotropy constant 
is proportional to J — K . At the isotropic point J — K 
the energy Wr = {K/2) J(yR)'^d'^x, and there exists a 
BP-type soliton that has the energy Ej^k ~ 2ttK and 
is a special case of the g = 1 CP^ solution p4|) . The 
CP^ charge q given by P^ is obviously equal to the 
Pontryagin index Qr defined by pT|) with m 1^ R; the 
BP solution corresponds to the mapping of S'^ onto the 
subspace CP^ embedded into CP^ and has Qr — q = 1. 
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FIG. 1: Energy of topological solitons in the vicinity of the 
S'?7(3)-symmetric point J — K. Solid lines show the energy 
of static solitons with the CP^ topological charge q — 2, and 
the dashed line corresponds to the descendants of the q = 1 
soliton of the CP^ model. 

For a finite "anisotropy" (J — K) the BP soliton be- 
comes unstable against collapse, but the situation is dif- 
ferent for the spin- nematic and FM regions. In the FM 
case (J > K) the anisotropy is of the easy-axis type, 
and there exist Qn = q = I dynamic solutions, with 
R precessing around the z axis 2], which are smoothly 
connected to the BP solitons in the J ^ K limit. A de- 
tailed analysis [9] shows that the minimal energy of such 
dynamic solitons exhibits a nonanalytical behavior of the 
type Emin = i^j=Ar(l + 3.74y^J/i^ - 1), as shown in Fig. 
[TJ For small J — K< O.IK, when the above expression is 
valid, the energy of a static q — 2 {Qm — 1) soliton con- 
sidered above stays higher than the energy of the q = 1 
[Qr — 1) dynamical soliton, but at the same time, it 
remains smaller that the energy of two dynamical q = 1 
solitons, which indicates that it is energetically favorable 
to bind two Qji — 1 solitons into a single Qm — 1 one. 

In the nematic case {J < K) we effectively have a 
ferromagnet with the easy-plane anisotropy. For such 
case, delocalized 7ri-solitons (vortices) exist. Vortices in 
i2-field correspond to spin-nematic disclinations consid- 
ered in Ref. [10]. The energy of a single vortex diverges 
logarithmically with the system size, so a static vortex- 
antivortex pair is unstable against collapse. The BP soli- 
ton can be considered as a pair of "merons" carrying 
topological charge = i each 11]. For small [K — J) 
those "merons" can be viewed as a vortex and antivor- 
tex with a finite out-of-plane component of the vector ii, 
they are subject to a gyroforce [10|, and there may ex- 
ist stable dynamic solutions (rotational pairs of vortices) 
similar to those studied in Ref. [H, [ij]. Their energy 
will tend to 2ttK in the limit J ^ K; similarly to the 
FM case, in the vicinity of the J — K point the Qr = 1 
topological solitons will be unstable against pairing into 
"nematic" Belavin-Polyakov solitons with Qu = 5- 

Finally, a few words are to be said about the other, an- 



tiferromagnetic 5'/7(3)-symmetric point J = 0. From (|4]) 
one can see that on any bipartite lattice the transforma- 
tion tj 1-^ t* for all j belonging to one sublattice maps the 
points J ~ and J = K onto each other. As can be seen 
from pH)) . and is especially clear from the "spin anal- 
ogy" (PO]) . for J < K/2 the short-range correlations are 
antiferromagnetic, and the proper transition to the con- 
tinuum description becomes more complicated; however, 
one can show that the difference concerns only dynamics 
and does not affect the static properties. The arguments 
leading to and thus the conclusion on soliton pairing 
equally apply to the vicinity of J = point. 

Summary.- We have studied the structure of topo- 
logically nontrivial solitons in a general non-Heisenberg 
model of the 2d isotropic 5 = 1 magnet. In the vicin- 
ity of special points with SU{S) symmetry the system 
can be described with the help of the CP^ model. It is 
shown that when the SU (3) symmetry is broken down to 
SU{2), solitons of the CP^ model with odd topological 
charge become unstable and bind into pairs. 
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